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ABSTRACT 

Four IRAS-detected early type galaxies were observed with the Infrared Space Ob- 
servatory. With the exception of the 15 fim image of NGC 1052, the mid-IR emission 
from NGC 1052, NGC 1155, NGC 5866 and NGC 6958 at 4.5, 7 and 15 /xm show ex- 
tended emission. Mid-IR emission from NGC 1052, NGC 1155, and NGC 6958 follows a 
de Vaucouleurs profile. The ratio of 15/7 /mi flux decreases with radius in these galaxies, 
approaching the values empirically observed for purely stellar systems. In NGC 5866, 
the 7 and 15 /im emission is concentrated in the edge-on dust lane. 

All the galaxies are detected in the [CII](158 /um) line, and the SOs NGC 1155 
and NGC 5866 are detected in the [01] (63 /mi) line as well. Previous detections of 
neutral ISM are sparse: only NGC 1052 had been detected in H I and NGC 5866 in 
CO. The ISO-LWS observations of the [CII] line are more sensitive measures of cool, 
neutral ISM than H I and CO by about a factor of 10-100. Comparison of [CII] with 
Ha shows that [CII] does not arise from HII regions and therefore must arise in the 
neutral medium. Three of four early type galaxies, namely NGC 1052, NGC 6958 
and NGC 5866, have low ratios of far- infrared to blue luminosity and show a lower 
L[ch]/Lfir, which is explained by postulating a softer radiation field from old stellar 
populations in early type galaxies, compared to spirals and irregulars, where young stars 
are present. While optical photons are effective in heating the dust, UV radiation is 
needed to heat the gas by the grain photoelectric mechanism. The low [CII] /CO ratio 
in NGC 5866 (Lj CII ]/Lco(l — 0) < 570) confirms this scenario. 

We estimate the UV radiation expected from the old stellar populations in these 
galaxies and compare it to that needed to heat the gas to account for the cooling 
observed [CII] and [OI] lines. In three out of four galaxies, NGC 1052, NGC 5866 and 
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NGC 6958, the predicted UV radiation falls short by a factor of 2-3 of that required to 
sufficiently heat the gas. In view of the observed intrinsic scatter in the "UV-upturn" in 
elliptical galaxies and its great sensitivity to age and metallicity effects, this difference 
is not significant. However, the much larger difference (about a factor of 20) between 
the UV radiation from old stars and that needed to produce the far-infrared lines for 
NGC 1155 is strong evidence for the presence of an additional UV source, probably 
young stars, in NGC 1155. 

1. Introduction 

A significant fraction of early type galaxies (ellipticals and SOs) contain some cool interstellar 
medium (ISM) including neutral gas and dust (Jura et al. 1987, Knapp et al. 1989). However, the 
physical conditions in the ISM in these galaxies are relatively unknown. Many ellipticals have been 
detected by the X-ray emission of hot gas (Roberts et al. 1991), and many in the H I and CO lines 
indicating cold-neutral medium (Lees et al. 1991, Wiklind, Combes & Henkle 1995), and through 
dust obscuration (Sadler & Gerhard 1985, Goudfrooij et al. 1994a 1994b). The environment in the 
ellipticals may lead to different equilibrium physical conditions in the various phases of the ISM 
because of the presence of hot X-ray gas and a higher ratio of stellar radiation compared to the 
cool gas and dust. But the IRAS colors in 60 and 100 micron emission indicate cold dust with 
temperatures comparable to found in spirals. 

We observed four elliptical/SO galaxies with the Infrared Space Observatory (Kessler et al. 
1996) to study the physical conditions in the gas and dust in these galaxies. These galaxies are 
part of a large sample of normal star- forming galaxies spanning a range of morphologies, infrared 
luminosities, FIR colors (F v (60)/F v (100)) and FIR-to-Blue ratio. ISO-CAM (Cesarsky et al. 1996) 
images at 7 and 15 microns are used to determine the distribution of mid-IR emission (section 3) 
and hence its origin. The mid-IR spectrum of NGC 5866 was taken with ISO-PHT (Lemke et al. 
2000) and is presented by Lu et al. (2000). The FIR spectroscopy of fine structure lines: [CII] (158 
fim) and [OI] (63 /im), was carried out with ISO-LWS (Clegg et al. 1996). Assuming [CII] and [OI] 
to be the major cooling lines, we estimate gas heating rates and examine possible heat sources in 
these galaxies in section 4. For NGC 1155 and NGC 5866 where both [CII] (158 jum) and [OI] (63 
fj,m) lines were detected we estimate the radiation and gas density, Go and n. 

The sample considered in this paper consists of two elliptical galaxies: NGC 6958 and NGC 1052 
and two SO galaxies: NGC 1155 and NGC 5866. Like any sample of four galaxies, these have varied 
properties. NGC 1052 harbors an active nucleus (Fosbury et al. 1978). NGC 5866 has a large edge- 
on disk. NGC 1155 is labelled an elliptical in the NASA Extragalactic Database (NED), but has 
been reclassified as an SO by Harold Corwin (private communication) based on morphology using 
B-band CCD images taken at Lowell Observatory. These images were taken with a CCD camera 
with 0.61" pixels and approximately 1.5" seeing. This galaxy is small (0.8 ' ) and distant which 



-3- 



may make a definitive classification difficult. The B-band surface brightness profile is well fit by a 
de Vaucouleurs profile and cannot be fit with a point source plus exponential disk. NGC 1155 is 
also labelled a starburst in NED and is a Markarian galaxy. It has a spiral companion (NGC 1154) 
two arcminutes away and a bridge can be seen between the two galaxies in deep optical images. 
NGC 6958 seems the most unremarkable galaxy in this set. 

2. Observations and data analysis 

The Mid-IR imaging of NGC 1155 and NGC 6958 was done with ISO-CAM at 6.75/«n (LW2 
filter, AA = 3.5/xm) and at 15 (LW3 filter, AA = 6/Ltm), using the raster scan mode to cover 
roughly 4.5' x 4.5' centered on the nucleus. CAM was set to 6"/pixel, and the raster was made 
up of 2 x 2 pointings separated by 81", or 13.5 pixels in each direction, allowing for better spatial 
sampling. At each pointing in the raster scan we integrated for 50 seconds in each band: 25 x 2 
second exposures at 15 /um, and 10 x 5 seconds at 7 fj,m. Typical sensitivity is 0.02 Jy/pixel in the 
7/ira filter and 0.04 Jy/pixel in the 15/xm filter. More details of CAM data reduction can be found 
in Dale et al (2000). NGC 1052 and NGC 5866 were observed by Vigroux et al. in the LW1 (4.5 
/xm), LW2 (7/im) and LW3 (15/zm) filters with 3 and 6 "/pixel. A detailed analysis of these data 
is presented by Madden et al. 2000. 

With the 80" beam of the long wavelength spectrograph (LWS) and spectral resolution of 
0.6 — 0.3//m we measure total line flux for these galaxies. The galaxies are estimated to have 
FWHM < 0.5' in FIR emission using deconvolved IRAS maps. The [CII] line observations were 
planned to achieve (la) sensitivities of 5 x 10~ 5 x Fpm, where Ffir is the total far- infrared flux 
of the galaxy and is computed according to the relation Ffir = 1.26 x 10~ 14 [2.58 x F l/ (60fim) + 
F u (100[im)]Wm~ 2 (Helou et al. 1988), where F„(60//m) and F^lOO/zm) are flux densities in Jansky 
at 60 and 100 fim measured by IRAS. 

The data were reduced and calibrated with the ISO data reduction pipeline OLP7.0. Post- 
pipeline data reduction and recalibration was carried out using an interactive data reduction package 
ISAP (cf. http://www.ipac.caltech.edu/iso/lws/lws.html). The line profiles were derived from 
several scans by running a median boxcar filter through the scans. We use the median of the 
observed fluxes instead of the mean to reduce the influence of outlying points arising from cosmic 
ray hits. Line fluxes were derived by integrating directly under the lines, after fitting a linear 
baseline to the continuum (cf Malhotra et al. 2000 for details on the LWS data reduction and line 
fluxes). 

Optical imaging of NGC 1155 was done at Lowell Observatory. The Ha images were obtained 
through narrow-band filters centered at the line; an off-band filter was used to image and subtract 
the continuum. Long slit spectroscopy and Ha and broad-band imaging of NGC 1155 was done 
at Lowell observatory. CO(1-0) observations on NGC 6958 were carried out at the Swedish ESO 
Submillimeter Telescope (SEST). 
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3. Mid-infrared imaging 

In spiral galaxies the mid-IR emission is dominated by small grains transiently heated to high 
temperatures and the fluorescence of large aromatic molecules (Puget & Leger 1989). This may 
not be the case for ellipticals given that the ISM-to-stellar ratio is low. A significant part of the 
emission at 12 microns could be from photospheres and circumstellar dust (Knapp, Gunn and 
Wynn- Williams 1992), although Sauvage and Thuan (1994) have argued that the IRAS color-color 
relation (Helou 1986) holds for early type as well as for late type galaxies indicating a similar 
origin from ISM. Hot dust near an active nucleus (e.g. AGN) can also be a major contributer in 
some cases. Because of its high spatial resolution and sensitivity, CAM imaging at 7 and 15 fim 
can tell us about the spatial distribution of the mid-IR emission in galaxies. The interpretation is 
complicated by the fact that the mid-IR emission maps depend on the distribution of the dust as 
well as the heating sources. 



Please see the accompanying 



JPG figure in the archive 



Fig. 1.— The mid-infrared and optical images of NGC 1052, NGC 1155, NGC 5866 and NGC 6958 
are presented for comparison. The images are (going left to right): DSS (digitised sky survey) 
images in R, 4.5 //m, 6.75/um and 15^m images from CAM. The 4.5/um images of NGC 1052 and 
NGC 5866 trace the starlight, while the 7/im and 15^m images of NGC 5866 trace the edge on dust 
lane. The 7/xm and 15/xm images of NGC 1052 are dominated by the active nucleus. All except 
the 15 fim image of NGC 1052 show extended emission. The orientation of the images is indicated 
by the arrows on the DSS images, where the arrowhead points north. The field of view is 4.5' on 
the side for NGC 1155 and NGC 6958 and 1.5' on the side for NGC 1052 and NGC 5866 

The mid-IR emission is extended in at least 9 out of 10 mid-IR images (Figure 1), although 
many of the images look compact because the light in ellipticals is very concentrated. Comparison 
of curves of growth (i.e. flux enclosed in a radius) show significant difference between a point source 
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and these galaxies (Figure 2) except for the 15 /urn image of NGC 1052. 

The 4.5 emission from NGC 1052 is extended and traces the stellar distribution. At 7 and 
15 /Ltm NGC 1052 comes closest to being dominated by the central point source, which is an active 
nucleus in this case. The flux in 7 and 15 micron bands in the central 6 " radius in NGC 1052 is 
about 80%, compared to 88% and 80 % expected for a point source in 7 and 15 micron filters with 
6 " pixel-scale. There is some uncertainty in the expected PSF depending on the placement of the 
point source within a 6" pixel. NGC 1052 was also observed in the 3 "/pixel mode on CAM which 
samples the PSF properly. The images with 3 "/pixel show that at 7/um NGC 1052 is extended 
and at 15 microns it is consistent with being unresolved (Figure 2). In this respect it is similar to 
NGC 3998 whose mid-infrared emission is dominated by the active nucleus (Knapp et al. 1996). 
The 7 and 15 (im. emission is extended for NGC 1155, where the central 6" accounts for only 47% 
and 56% of the total flux at 7 and 15 fim. In NGC 5866 the central 6" contains 24% and 40% of 
the total flux in 7 and 15 /um bands. In NGC 5866 the mid-IR emission at 7 and 15 fim traces the 
nearly edge-on dust lane, whereas the 4.5/xm emission traces the stars (Figure 1). 

To increase the signal-to-noise in the outer parts of the galaxies we average the emission in 
annuli after subtracting the background determined at radii greater than 90". Figure 3 shows the 
radial profile of the surface brightness in three galaxies, NGC 1052, NGC 1155 and NGC 6958. 
The radial surface density profile is derived by averaging azimuthally. The error-bars in surface 
brightness in Figure 3 are derived from the scatter in pixel brightness in a given annulus. A 
de Vaucouleurs radial profile provides a reasonable fit to the mid-IR surface brightness profile. The 
effective radius for the mid-IR emission is intermediate to that seen in optical B-band and Ha light 
for NGC 1155 (Figure 4). We also note that in all three galaxies the 15^m/7^m ratio decreases 
going from the center of the galaxy outwards (Figure 3). 

The 15^m/7^m ratio can be used as a further diagnostic to determine the origin of the mid- 
infrared emission in the current sample. In quiescent spiral galaxies this ratio is 1 (e.g. NGC 6946, 
Helou et al. 1996). In more actively star-forming regions/galaxies the 15/(/m/7^m ratio increases 
due to higher heating radiation (Vigoroux et al. 1999, Dale et al. 1999). 15^m/7/(/m < 1 is 
found in early type galaxies where the 7fim band contains a significant contribution from stellar 
photospheres (Madden et al. 1997). In the inner parts of NGC 1052 and NGC 1155, 15^m/7/xm is 
a factor ~ 2 higher than the ISM in quiescent spiral galaxies (Figure 3). This is presumably due to 
higher radiation density in the central regions of these galaxies due to higher density of stars or the 
active nucleus in the case of NGC 1052. The decrease in 15^m/7//m in the outer parts of the galaxy 
is due to different proportions of photospheric and ISM contributions, with the photospheres having 
a lower 15^m/7^m ratio (i.e., bluer emission). This scenario is consistent with the 15;um/7^m ratio 
of 0.5 seen in the outer parts of NGC 6958. This ratio is within the range of observed ratios in 
ellipticals devoid of ISM (Madden et al 1997), and indeed NGC 6958 does not show extended Ha 
emission beyond its nucleus, indicating that the ISM is not very extended in this galaxy (Phillips 
et al. 1986) . 
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The spatial distribution of the mid-IR emission could also reflect the distribution of sources 
heating the dust. In case of NGC 1155 the mid-IR emission follows a de Vaucouleurs profile with an 
effective radius R e which is intermediate between Ha and B band indicating that both ionizing and 
non-ionizing radiation heat the dust. A similar result was found for the spiral galaxy NGC 6946 
(Malhotra et al. 1996). If so, we would predict that the dust is quite extensive in NGC 1155. While 
the extent of the ISM is not known for NGC 1155, the luminosity in the [CII] line is as much as a 
third of the Milky Way value indicating the presence of a significant amount of gas. The presence 
of an extended component of dust can explain the discrepancy between dust masses obtained from 
optical extinction and FIR emission, the former being sensitive only to dust concentrated near the 
center (Goudfrooij et al 1994). 



4. Far-Infrared Spectroscopy 

The physical conditions in the different phases of the ISM in early type galaxies are not very 
well known, partly because there is not as much cool ISM in these galaxies, and partly because we 
do not live in such a galaxy. Much of our detailed knowledge about the various phases of the ISM 
in spirals comes from being in these phases and able to do absorption studies. In ellipticals/SOs 
hot gas is detected in the X-ray, and ionized gas is detected by Ha emission. The cold neutral 
medium is detectable as H I (21 cm) and CO by present instrumentation for typical masses of 
10 7 — 10 s Mq. As a result the detection rate for cold gas is about 15-50% (van Gorkom 1997, Rupen 
1997) depending on the morphological classification of galaxies and preselection of samples. The 
detection rate is higher for Ha (about 50%) where one is sensitive to 1O 5 M of gas (Phillips et al. 
1986). 

The far infrared fine structure lines dominate the global cooling of the neutral ISM of galaxies. 
Of these, [CII] (158 fiia.) and [01] (63 /jm) are generally the strongest, with [CII] dominating the 
lower densities and moderate UV radiation field intensities (Hollenbach, Tielens & Takahashi 1991; 
Hollenbach & Tielens 1997, 1999). We detect the [CII] transition in all four galaxies, albeit at low 
(3-<j) S/N for NGC 6958. The other major fine structure cooling line, [01] (63 /um), is detected 
in NGC 1155 and NGC 5866. In Table 1 we summarize what is known about the warm and cool 
ISM in the four early type galaxies. Detections are few - H I (21 cm) is detected in NGC 1052 
(van Gorkom et al. 1986), CO in NGC 5866 (Thronson et al. 1989) and Ha is detected in all four. 
There is tentative detection of CO in absorption against the nucleus in NGC 1052 (Knapp & Rupen 
1996). In Table 1 Mh(C + ) is the minimum mass of the hydrogen associated with the observed 
C+ emission and is derived in the high density, high temperature limit (n# > 3000, T 3> 91-fT), 
assuming solar abundance of carbon with all the carbon in the form of C+ (Crawford et al. 1985). 
Mh{C + ) is typically 10-100 times smaller than the observed detection or upper limits on HI and 
H 2 . 
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Table 1: Mass in Interstellar components 



Name 


Morphology J 


Distance 


M(HI) 


M(H 2 ) 


M(HII) 


M H (C+) 


Mdust 1 








(Mass 


in 








NGC 1052 


E4 


14 Mpc 


0.7 x 10 7 a 


< 3 x 10 7 b 


3 x 10 4 c 


> 4 x 10 5 


4.8 x 10 4 


NGC 1155 


SO 


45 Mpc 


? 


? 


? 


> 1 x 10 7 


1.8 x 10 6 


NGC 5866 


SO 


11.4 Mpc 


< 3 x 10 7 f 


> 3 x 10 8 s 


4 x 10 4 h 


> 9.7 x 10 5 


1.2 x 10 6 


NGC 6958 


E+ 


28 Mpc 


< 3 x 10 s d 


< 9 x 10 8 c 


? 


> 8 x 10 5 


3.5 x 10 5 



"from van Gorkom et al. 1986, adjusted to the LWS beam-size of 70" 

6 Knapp & Rupen 1996, Wiklind et al. 1995 

c Plana & Boulesteix 1996 

d Walsh et al. 1990 

e recent observations at SEST 

•fHaynes et al. 1990 

9 Thronson et al. 1989 

h Plana et al. 1998 

'Roberts et al. 1991, adjusted for the given distances 

J from NASA Extragalactic Database and Harold Corwin (private communication) 

4.1. Trends in line strengths in Early type galaxies 

The [CII] line strengths of normal spiral and irregular galaxies follow well defined trends 
with far-infrared colors (F iy (60 / um)/F^(100 / um)) indicating dependence on dust temperature and 
therefore on dust heating radiation density; and on star-formation activity, measured by the ratio 
of far-infrared and blue band luminosity Lfir/Lb- The ratio L^hj/Lfir decreases dramatically 
with increasing F^(60^m)/F^(100 / um) and L Fm /L B (Malhotra et al. 1997, 2000). 

While [CII] and [OI] are the main coolants of warm (~ 90 K), neutral gas; the main source of 
heating are the photoejected electrons from dust grains (Watson 1972). Only UV light is effective 
in producing photoejected electrons since the potential barrier for neutral grains is approximately 6 
eV, or more if the grains are positively charged (Bakes & Tielens 1994). The far-infrared continuum 
reflects the cooling, and thus the total heating of the dust grains. For galaxies with the most active 
star-formation and warm dust and hence most intense radiation fields, the dust grains become 
positively charged making the photoelectric heating less efficient. Therefore we see a decrease in 
[C II] and [O I] lines strengths relative to the dust continuum for the most active galaxies with 
higher F^(60/xm)/F^(100/xm) and Lfir/Lb (Figure 5). 

On the other extreme, galaxies with low rates of current star-formation, and softer radiation 
field a low Lrcnj/LpiR is expected since only the UV photons (A < 2000^4) contribute to photoelec- 
tric heating of gas, while dust is heated by optical as well as UV photons. In early type galaxies, 
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which typically have lower active star-formation rate relative to the total blue luminosity, we do 



In Figures 5a and 5b we plot these trends for the whole sample of distant galaxies with the early 
type galaxies circled. It is apparent in these figures that the four Es and SOs lie within the broad 
correlation between Lpuj/LpiR and F„(60//m)/F„(100//m) and Lfir/Lb. NGC 6958 and NGC 1052 
and NGC 5866 lie at the low end of Lfir/Lb, indicating a low rate of current star-formation and 



The [CII] emission could be from ionized gas in dense HII regions, diffuse ionized gas, neutral 
atomic medium, or from UV irradiated surfaces of molecular clouds. In nuclear regions of galaxies 
[CII] arises mostly from dense PDRs (Stacey et al. 1991, Crawford et al. 1985). However, integrated 
over the disks of normal spiral galaxies, a significant fraction (up to 50%) may also arise from diffuse 
medium, i.e. the cold neutral medium (CNM) (Madden et al. 1993) or from extended, low density 
ionized gas (Heiles 1994). With the limited knowledge of the phases of the ISM in these galaxies and 
without assuming that the temperature and density in these phases is very similar to the Galactic 
values, it is hard to identify the phase of the ISM which produces most of the [CII] emission. 

We can rule out HII regions as sources of much [CII]. Assuming case B recombination and 
T e = 10 4 K, one can place an upper limit on [CII] emission from HII regions (Malhotra et al. 2000) 
as 



here xc is the gas phase abundance of carbon with respect to hydrogen, and n e is the density of 
electrons in the HII region in cm -3 . The inequality sign is due to the fact that not all of the C in 
HII regions exists as C + ; some will be in the form of C 2+ . 

For low electron density (n e — ► 0): 



see a lower L[ CII ]/L F i R . 



low extinction. The Ljchj/Lfir ra ti° turns down slightly at the low end of star-formation. 



4.2. Origin of [CII] 




Lcii 
Lh« 



< 0.33x_4 



where x_4 = xc/W 4 In the high density limit (n e >> 17cm 3 ): 




For NGC 1155, NGC 5866 and NGC 6958 the ratio L([C 1 1]) / L(H a) is measured to be 2.5, 
0.5 and 58 respectively indicating that most of the [CII] is not associated with HII regions. For 
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NGC 1052 the observed L([C 1 1]) / L(H a) = 0.16. Taking n e = 10 3 cm~ 3 , inferred from the ra- 
tio of SII lines (Kim 1989), the contribution of HII regions to CII emission is expected to be 
L([CII])/L(Ha) < 0.015 (for x_4 = 3); an order of magnitude lower than the observed levels. 
Thus it appears that not much of [CII] emission arises from the ionized gas in early type galaxies. 
The total extinction is low so the Ha flux is unlikely to be extincted by more than a factor of two. 

Could the [CII] emission be coming from diffuse HI? It is hard to answer this question given that 
HI column density has been measured only in NGC 1052. For a cold neutral medium (CNM) with 
typical Galactic values T=80K and n H = 90cm~ 3 , M{HI)/L{CII) = 20(3 x W~ 4 /X C + )M Q /L Q . 
For NGC 1052 M(HI)/L(CII) = 14 which can easily be accommodated with a slightly higher gas 
density or warmer gas, or some contribution to [CII] from molecular medium. It is quite likely that 
most of the [CII] emission is associated with molecular gas in NGC 5866, given that the upper limits 
on H I mass are ten times lower than the observed M(H2). Unfortunately, no useful constraints 
can be placed on the origin of [CII] in NGC 1155 and NGC 6958 where neither H I nor CO have 
been detected. It should be easy to estimate the H I contribution to the origin of [CII] emission in 
NGC 1155 with H I observations. If [CII] is associated with H I we should see M(HI) ~ 1O 8 M . 



4.3. [CII] /CO ratio 

For [CII] emission arising from photo-dissociation regions, softer radiation leads to a lower 
fraction of gas in C + due to reduced photodissociation of CO molecules, making a thinner C + 
layer on the surface of the molecular cloud (Spaans et al. 1994). This scenario predicts a lower 
ratio of the flux in [CII] versus various CO transitions as well as lower L[ch]/Lfir. For Galactic 
HII and OB star-forming regions and for starburst galaxies the ratio L[cn]/Lco(l — 0) = 4100, 
whereas for more quiescent spiral galaxies this ratio was found to be L[cii]/Lco(l — 0) = 1500 
(Stacey et al 1991). For NGC 5866, which is the only galaxy in the current sample with a CO 
measurement Lpn]/Lco(l — 0) < 570. This is an upper limit because the CO is measured in 
the 45" beam of FCRAO (Thronson et al. 1989) and the [CII] is measured in a 70" beam. This 
provides evidence for a soft radiation field. For NGC 6958 and NGC 1052 we only have upper 
limits on CO(1-0) emission. For NGC 1052 L[ CII ]/Lco(l — 0) > 1100, which is comparable to 
the quiescent spiral galaxy sample, but most of the [CII] seems to be associated with the atomic 
gas in this galaxy, while the low Lrcni/Lco(l ~~ 0) prediction in soft radiation fields is made for 
dense PDRs on the surfaces on molecular clouds. In NGC 5866, where we find low Lrcn]/Lco(l — 0) 
values, the M(HI)/M(H2)< 0.1 indicating that most of the [CII] is associated with molecular clouds, 
providing a cleaner diagnostic of the softer radiation field. For NGC 6958 upper limits on CO imply 
L[cii]/Lco(l — 0) > 170 and more sensitive CO observations are needed. 
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4.4. Physical conditions in the gas 

With the few FIR line measurements of these galaxies we can attempt to derive the physical 
conditions in the warm, neutral, FIR emitting gas. 

NGC 1155: We compare the two FIR line ([CII] and [OI](63//m)) and line to continuum 
ratios in NGC 1155 to models of photodissociation regions by Kaufman et al. 1999. Using the 
ratios ( [CII] + [01]) /FIR and [OI]/[CII] gives density n ~ 10 2 cm- 3 and FUV flux G ~ 10 2 (G 
conventionally is the FUV flux normalized to the average local interstellar flux of 1.6 x 1CT 3 ergs 
cm -2 s _1 (Habing 1968)) There are no HI or CO measurements for this galaxy. The chances of 
detecting both are excellent, both because of the high minimum mass of hydrogen associated with 
the C + emission and because of the density of the gas derived. 

NGC 1052: The column density of HI is measured to be 2.5 x 10 20 cm~ 2 (van Gorkom et 
al. 1986) near the nucleus and about 1 arcminute away. It is then reasonable to assume that 
N(HI) is constant in the LWS 70" beam and the [CII] emitting gas fills the beam: I([CII])= 
1.1 x 1CT 6 ergs / cm 2 / s / sr . Using this in the following equation with the further assumptions that 
the temperature of the gas T > 91K and the ionization fraction is small (X e < 1CT 3 , so we can 
neglect collisional excitation due to electrons), we can solve for the density of the gas njj. 

_2i ~ , 2e(- 91 /r) 



I {[CI I}) = 2.3 x 1CT 21 N{HI) X c + ( 



l + 2e(-^m +n crit /n H ' 



Where n cr it = 2 x 10 3 is the critical density for collisions with HI (Launay & Roeff 1977, Hayes 
& Nussbaumer 1984) and X^+, the C + abundance is taken to be the Galactic value of of 3 x 10 -4 . 
The density of the HI then is derived to be njf = 6cm~ 3 , consistent with [CII] arising from diffuse 
HI. 

NGC 5866 Comparing the models to the observed ( [CII] + [01]) /FIR, [CII]/[OI] and [CII]/CO(l- 
0) ratios does not yield a solution for gas and radiation density in NGC 5866, because the models 
assume a harder radiation field, and that half the dust heating is from far-UV. In softer radiation 
fields from old stars the optical light contributes a larger fraction of the dust heating and therefore 
to the FIR continuum. 

Ignoring ( [CII] + [01]) /FIR ratio, we estimate the intensity of [CII] emission in the beam, as- 
suming a beam filling factor of 1, to be I([CII])= 5.5 x W~ 6 ergs/cm 2 /s/sr. Using this and the ratios 
[01] /[CII] =0.2 and [CII]/CO(1-0)< 570, PDR models of Kaufman et al. (1999) give the density 
n = 3 x 10 3 cm~ 3 and UV radiation density Go = 1 — 3. The same models estimate a temperature of 
about 30-40 K at the surface of the molecular clouds and a thermal pressure of ~ 10 5 i"Ccm~ 3 , which 
is 30 times higher than the local solar neighborhood value, but 10 times lower than the pressure at 
the center regions of the Milky Way (10 6 Kcm~ 3 : Spergel & Blitz 1992). This is a plausible, since 
the LWS beam encompasses gas at radius < 2 kpc from the center of NGC 5866. The PDR models 
then overpredict ( [CII] + [OI]) /FIR = 5 x 10~ 3 , when the observed ( [CII] + [OI]) /FIR = 2 x 10~ 3 . 
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The PDR models assume that equal amounts of heating of the dust comes from UV and non-UV 
photons. The current discrepancy in ([CII] + [OI])/FIR ratio thus implies that four times as much 
heating of the dust grains comes from non-UV radiation as from UV. This is also our favored 



In spiral galaxies [CII] is frequently used as a measure of star-formation activity (Stacey et al. 
1991). In early type galaxies emission in the fine structure lines [CII] and [01] do not necessarily 
indicate current star-formation (or PDR emission). In spirals, [CII] and [01] are the major coolant 
of warm, neutral gas and most of the heating is by photoelectrons from dust grains (Watson 1972). 
The potential barrier to photoejection from neutral grains is approximately 6 eV (Bakes & Tielens 
1994), or higher if the grains are positively charged. This means that only photons shortward of 
~ 20004 contribute to heating of the gas. Old stellar populations (PNe, HB, PAGB stars) in 
Es/SOs produce enough ionizing flux shortward of 912 A to explain the Ha emission seen (Rose & 
Tinsley 1974, Binette et al. 1994) 9 . So it is possible that they contribute to heating the gas by 
photoelectric heating. But what is the magnitude of that contribution? 

To answer that question we estimate photoelectric gas heating by UV photons from old stellar 
populations. We use GISSEL (Bruzual and Chariot 1993) stellar population models of instanta- 
neous burst of ages 10 and 13 Gyr to calculate the UV spectra of early type galaxies. The mass 
of each galaxy is normalized to match the B-band flux within the LWS beam of 70". An upper 
limit to the gas heating rate is calculated by assuming that the work function of the grains is 6 eV 
and the yield of photoelectrons heating process is 7 = 0.1 (Tielens & Hollenbach 1985). This is an 
upper limit on the heating and a lower limit on the UV required because it neglects the positive 
charging of the grains by photoelectric process, which raises the potential barrier for ejecting elec- 
trons off grains, therefore requiring even more energetic photons. We neglect here the contribution 
of negatively charged grains, which have an ionization potential as low as 2 eV, because they are 
estimated to contribute less than 10% of the heating (Spaans et al. 1994). The total power that 
goes into heating the gas is given by: 



In the above equation N v are the number of photons at frequency v and f a b s is the fraction of 
photons absorbed by the dust. f a t, s should be less than unity given the low ratio of ISM to stars. 
fabs is estimated from the observed values of Lb /Lfir and using a standard extinction curve with 



scenario for the observed low L[ CII ]/Lfir and L[ CI r]/Lco(l — 0). 



4.5. Gas energetics 




9 If the ionizing photons come from old stellar populations, one might expect the Ha distribution to be diffusely 
distributed, but it is usually seen to be associated with dust lanes (Goudfrooij et al. 1994) 
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the spectral energy distribution from stellar population models. We note here in passing that it was 
not possible to reproduce the high Lfir/Lb = 0.29 seen in NGC 1155 with any of the old population 
models by any extinction values. For the other three galaxies Ay varied from 0.05 to 0.25. Table 
2 lists the heating luminosity Lpe from old stellar populations, compared to the luminosity of the 
[CII] and [01] lines; for NGC 1052 and NGC 6958 the [OI] line strength is estimated using the 
correlation between [OI]/[CII] and F I/ (60/xm)/F I/ (100/im) (Malhotra et al., 2000). 

The UV heating produced by the old stellar populations fails to account for the [CII] and [OI] 
emission by factors of 2 or 3 in NGC 1052, NGC 6958 and NGC 5866, and clearly fails for NGC 1155 
by a factor of about 20. The theoretical estimates of UV production by old populations are uncertain 
and vary significantly with age and metallicity (cf. Magris and Bruzual 1993). Therefore a factor of 
2-3 discrepancy in the required UV and theoretical UV production does not necessarily imply the 
presence of extra UV sources. However, for NGC 1155 we can conclude that the strength of [CII] 
and [OI] emission imply a higher UV production than can be explained by the old populations. 

This conclusion is also supported by the optical spectrum of NGC 1155, which shows narrow 
emission lines. The ratio of [Nil] (6583 A)/Ha ~ 0.5 characteristic of HII regions (Figure 6). 
NGC 5866 shows a ratio of [Nil] (6583A)/Ha ~ 1.34 and Ho et al. (1997) classify the nucleus as a 
transition object between a LINER and an HII nucleus. 

The observed far-UV-optical colors (1550-V) of Es/SOs vary by a factor of 10 for active and 
quiescent galaxies (Brown et al. 1998, Burstein et al. 1988, hereafter B3FL). The scatter is reduced 
by using a correlation with the magnesium index Mg2. Quiescent E/S0 galaxies, with or without 
an active nucleus, show a higher Mg2 index than galaxies with star-formation, except when they 
have very red (1550-V) colors. We compare our sample to B3FL figure 1(a), except we estimate the 
flux in the 1500A band from the UV required to produce the [CII] and [OI] line emission (hereafter 
referred to as FUVp//]), since 3 out of 4 galaxies (except NGC 1052) do not have their far-UV 
measured. We cannot comment on NGC 1155 whose V or Mg2 measurements are not available. 
NGC 6958 and NGC 5866 have low Mg 2 indices, 0.23 and 0.21 respectively (Davies et al. 1987, 
Huchra et al. 1996, Fisher et al. 1996). NGC 5866 has normal FUV [ar/] -V color so the low Mg 2 
could be due to recent star-formation, or due to extinction by the central dust lane (Fisher, Franx 
& Illingworth 1996). NGC 6958 has red color FUVp//]-V=4.2, and lies on the locus of quiescent 
galaxies with red colors and low Mg 2 indices, indicating little recent star-formation. For NGC 1052 
the FUV[ C 77]-V=1.08 magnitudes, 2 magnitudes bluer than the IUE observed color (B3FL). In 
other words, from the [CII] estimate NGC 1052 is bluer than the sample of quiescent galaxies. The 
FUVr^jn is also 6 times higher than the observed UV in this galaxy, although some uncertainty 
remains because IUE measured fluxes between 1300 and 3000 A and FUV[c//] is an estimate of flux 
between 912 and 2070 A . The extra UV in NGC 1052 could come from very obscured source(s), 
either star-formation or the central AGN. 

Other sources of gas heating include X-rays from hot gas and active nuclei. Only NGC 1052 
shows evidence for a Seyfert nucleus and is detected in X-rays by Einstein and ROSAT. It shows a 
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higher luminosity in X-rays than the average Lx~Lb relation (Fabbiano et al. 1992) and the X-ray 
gas temperature is much higher than expected from the average relation between stellar velocity 
dispersion and temperature of the X-ray gas, indicating that X-ray emission is dominated by AGN 
or X-ray binaries rather than diffuse hot gas (Davis and White 1996). The X-ray luminosity of 
NGC 1052 is ~ 3 x 10 6 L Q . About 10% of X-ray luminosity could be used for gas heating in 
atomic medium and about 40% in molecular medium in the so called X-ray dissociation regions 
XDRs (Maloney et al. 1996). In NGC 1052 X-rays could make a significant contribution to gas 
heating. But in XDRs one expects a higher line to continuum ratio for the [CII] and [01] lines 
because the X-rays are equally efficient at heating the gas and the dust, whereas the UV radiation 
in PDRs heats dust much more efficiently Because NGC 1052 shows a lower L[ CII ]/Lfir than the 
average for other (mostly spiral) galaxies in the sample, a significant X-ray heating of the cool 
gas seems unlikely. NGC 1155 and NGC 6958 were not detected by Einstein or ROSAT all sky 
survey. We estimate the X-ray luminosity from the Lx — Lb relation (Fabbiano et al. 1992) to 
be ~ 1O 5 L and ~ 0.7 x 1O 5 L for NGC 6958 and NGC 1155 respectively. That is simply not 
sufficient to heat the gas to produce the [CII] flux seen (Table 2). NGC 5866 has a luminosity of 
Lx(A^5866) = 0.7 x 1O 5 L , again insufficient to heat the cool gas enough to explain the flux in the 
far- infrared cooling lines. 



5. Summary and Conclusions 

From observations of two elliptical galaxies and two SOs with ISO-CAM and ISO-LWS we can 
conclude that 

(1) The Mid-infrared emission from early type galaxies arises both from ISM and photospheres 
of stars. 7/15 ftm ratio varies with radius indicating that ISM dominates at the centers and in 
NGC 6958 is not very extended. In NGC 1155 the dust emission extends to 10 Kpc. In NGC 1052 
the Seyfert nucleus dominates the mid-IR emission. In NGC 5866, the mid-IR emission at 7 and 
15/im follows the central dust lane. 

(2) [CII] (158 fim) and [OI] (63 fxm) lines are detected in galaxies where H I and CO detections 
were difficult, making ISO-LWS one of the sensitive probes of small quantities of cool ISM (T ~ 100 
K). [CII] is not mostly from classical HII regions. [CII] emission is associated with H I in NGC 1052, 
and with molecular gas in NGC 5866. 

(3) The line to continuum ratio L[ CII j/Lfir measures the efficiency of gas heating by pho- 
toelectric heating. In three of the four early type galaxies considered here L[ CII ]/LpiR is lower 
by a factor of 2-5 than the typical values in a sample of 60 normal galaxies. A softer radiation 
field, i.e. relatively deficient in UV is the most likely explanation. This is corroborated by a low 
L[cn]/Lco(l — 0) ratio observed in NGC 5866. 

(4) The fluxes in the cooling lines, [CII] and [OI], provide fairly stringent lower limits on the 
UV flux between 912 and 2100A With all the uncertainties in modelling, UV from evolved stellar 
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populations is insufficient by a factor of 2-3 for NGC 1052, NGC 5866 and NGC 6958 and a clear 
factor of 20 for for NGC 1155. Heating of cool ISM by X-rays is a plausible scenario only for 
NGC 1052. In NGC 6958, NGC 1052 and NGC 5866 it is plausible that old stars produce enough 
UV to heat the ISM to explain the cooling by [CII] line. In NGC 1155, young stars are needed to 
provide the UV heating to account for the cooling via the FIR lines. 

(5) Comparing the [CII] and [01] (63/xm) line flux measurements of NGC 1155 with PDR 
models by Kaufman et al. 1999 we infer gas density n ~ 10 2 cm~ 3 and UV radiation density 
~ 10 2 Go, where Go is the UV radiation density in the solar neighborhood. Comparing the 
models to the observed [CII], [OI] and CO(l-O) line ratios for NGC 5866 yields G = 1 - 3, 
and n = 3 x 10 3 cm -3 , and the estimate that non-UV radiation contributes about 80% of the dust 
heating in this galaxy. 
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Fig. 2. — This plot shows the cumulative flux profiles of elliptical galaxies at 7 and 15/xm as 
compared to the point source profiles. NGC 1155 and NGC 6958 were observed in a mode with 
1 pixel=6 " and they are plotted in the lower two panels. NGC 1052 was also observed in the 3 
" pixel mode and its profile is plotted in the upper two panels. Except for the 15 fim image of 
NGC 1052, all other images are clearly extended. This is manifest from the slow rise of cumulative 
flux with radius for galaxies as compared to the point source profile. Within measurement errors 
the 15 fim profile of NGC 1052 is consistent with an unresolved source. 
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Fig. 3. — The surface brightness of early type galaxies in the mid-infrared at 7 and 15 microns 
follows de Vaucouleurs radial profiles (left and middle columns). The ratio of LW3 (15 [im) and 
LW2 (7 /im) light decreases towards the outer parts of galaxies (right column). The 15/um/7^m 
ratio for NGC 6958 reaches an asymptotic value of 0.5, characteristic of old stellar populations in 
galaxies devoid of ISM. This indicates that the ISM is concentrated near the center. 
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Fig. 4. — The radial profile of optical B band emission (filled circles), the Ha emission (open circles), 
and the mid-infrared emission (7/xm: open squares, 15/xm: triangles) in NGC 1155 is plotted along 
with the best fit de Vaucouleurs profile for each wavelength. The mid-infrared light shows a radial 
distribution intermediate to the B band and the Ha. 



Table 2: Luminosities in cooling lines compared to photoelectric heating 
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NGC 1155 
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NGC 6958 
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(0.66) b 
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0.46 
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a Fluxes in the [OI] (63 fim) and [CII] lines are taken from Malhotra et al. 2000, typical uncertainty in a line flux is 
30% 

6 [OI] (63 /im) flux is derived from the [CII] flux and the IRAS colors F v (60) / F v (100) using the correlation seen between 
[OI]/[CII] and F I/ (60)/F„(100) (Malhotra et al. 2000). 
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Fig. 5. — The ratio of FIR and [CII] line luminosities is plotted against the dust temperature 
characterized by the ratio of flux densities at 60^m and 100/mi from IRAS measurements for a 
sample of galaxies spanning a whole range of morphologies (cf Malhotra et al. 2000). The four 
early type galaxies are shown with the bull's eye symbol, and lie within the broad trend of decreasing 
L[ CII i/Lfir with the increase in F iy (60/um)/Fj y (100 / um). (b) L[ CII ]/Lfir shows a similar trend with 
the ratio of FIR/Blue luminosity, which can be used as an indicator of star formation activity. Three 
out of four early type galaxies are at the lower end of mass normalised star-formation activity. These 
galaxies-NGC 1052, NGC 6958 and NGC 5866 show a lower L[ CII j/Lfir due to softer radiation field 
from aging stellar populations. 
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Fig. 6. — The optical spectrum on NGC 1155 shows narrow lines and line ratios characteristic of 
HII regions; i.e. Ha/[NII] >1) lending support to our conclusion from the analysis of FIR lines 
that there has been recent star-formation. 
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